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Abstract 


We use numerical and analytical tools to demonstrate arguments in favor of the existence of 
a family of smooth, symplectic diffeomorphisms of the two-dimensional torus that have both a 
positive measure set with positive Lyapunov exponent and a positive measure set with zero Lya¬ 
punov exponent. The family we study is the unfolding of an almost-hyperbolic diffeomorphism 
on the boundary of the set of Anosov diffeomorphisms, proposed by Lewowicz. 

1 Introduction 

Due to an extremely complicated intermixture of regular and chaotic orbits, the problem of the orbit 
structure of a generic, smooth symplectic map remains mainly open, even for the two-dimensional 
case. When the map is sufficiently smooth, its phase space typically exhibits both regular dynamics 
due to invariant KAM curves (for instance, in the neighborhood of elliptic periodic orbits) and seas 
of chaotic orbits (which numerical investigations indicate can be densely covered by a single orbit). 
Moreover, such structures are observed—again in numerical simulations—to occur at all scales. All 
this is well known and shown in many papers, for a review see, e.g., [Mei92]. It is generally agreed 
that no tools currently exist that allow one to rigorously elucidate the main points of this observed 
picture [Sin95]. Of course, selected parts of this landscape can be explained; for example, KAM 
theory provides a proof of the existence of invariant curves near generic elliptic periodic points. 
However even for this case, there is essentially no rigorous characterization of the orbit behavior in 
the so-called chaotic zones—as depicted in Arnold’s famous sketch [Arn63, Ber78]. 

*LML was supported in part by RFBR (grant 14-01-00344a), the Russian Ministry of Science and Education 
(project 1.1410.2014/K, target part) and the Russian Science Foundation (project 14-41-00044). Useful conversations 
with V. Grines are gratefully acknowledged. Both authors acknowledge a support from CRDF (grant FSAX-14-60273- 
0 ). 
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There has been much study of the destruction of invariant curves, and the resulting transition 
from regular (quasiperiodic) to irregular (chaotic) behavior, in parameterized families of area¬ 
preserving maps. Since a smooth invariant curve is not isolated, its destruction is caused by a 
loss of smoothness and, at least for twist maps, to the formation of a new, quasiperiodic invariant 
Cantor set; an Aubry-Mather set [ALD83, Mat82]. In many families, one observes the ultimate 
destruction of all the invariant circles (of a given homotopy class), and this leads to the study of 
the “last” invariant curve, and the development of Greene’s residue criterion and renormalization 
theory [Mac93]. 

At the opposite extreme, the ergodicity and hyperbolicity properties of Anosov diffeomorphisms 
are well-understood [Fra70]. This extreme of uniformly hyperbolicity can be thought of as a com¬ 
plementary limit to integrability: the study of perturbations from “anti- integrability” was initiated 
in [AA90]. Aubry’s results are based on the consideration of infinitely-degenerate diffeomorphisms 
and provide proofs of the existence of horseshoes; however, they do not lead to proofs of a positive 
measure of chaotic orbits. 

There have been attempts to understand the dynamics of symplectic diffeomorphisms on the 
torus beyond the boundary of the Anosov maps [Prz82]. Przytycki proved the existence of a 
curve of diffeomorphisms that cross the Anosov boundary such that, outside the boundary, there 
is a domain on the torus bounded by a heteroclinic cycle formed by merged separatrices of two 
saddles that contains a generic elliptic fixed point. The remaining set of positive measure has 
a nonhyperbolic structure and positive Lyapunov exponent. The drawback of this example is in 
its infinite codimension in the space of smooth symplectic diffeomorphisms with C^-topology: the 
merging of separatrices of saddles is a codimension-infinity phenomenon. Przytycki’s family unfolds 
a smooth, almost-hyperbolic symplectic diffeomorphism of the torus proposed earlier by Lewowicz 
[Lew80]. This diffeomorphism is a K-system that has positive Lyapunov exponent [CEOl]. 

This same trick (with the same drawbacks) was used later in [Liv04] to construct symplectic 
diffeomorphisms arbitrarily close to Lewowicz’s almost-hyperbolic map. Smooth symplectic, transi¬ 
tive diffeomorphisms that are K-systems on closed two-dimensional manifolds other than tori were 
constructed in [Kat79] (see also, [GK82]). Again it is not clear how these results can be used to 
understand the orbit structure for a generic diffeomorphism. 

Following [LerlO], we study the map / : —?• T^, where T = M/Z, defined through 

f{x, y) = {x + y + g{x),y + g{x)) mod 1. (1) 

If the “force” g were a degree-zero circle map, then (1) would be a generalized Chirikov standard 
map [Mei92]. Instead, we assume that g is a degree-one, circle map: 

g{x + 1) = g{x) + 1. (2) 

When 5 is a monotone increasing diffeomorphism, (I) is Anosov: every orbit is uniformly hyperbolic 
and / is topologically conjugate to Arnold’s cat map a : —)■ T^, 

/2 l\ 

mod 1, where ^ ~ ^ ^ 

More generally, Franks showed that (1) with (2) is semi-conjugate to a [Fra70], i.e., there is a 
continuous, onto map A: : —>■ such that 

k o f = a o k. (4) 


a{x,y) = A 
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The map k depends continuously on g, and when g is strictly monotone, A: is a homeomorphism, 
implying—as mentioned above—that / is then conjugate a. 

In [LerlO], the first author made an attempt to elucidate the features of (1) when the circle 
map g acquires a critical fixed point, 


g{xp) = Dg{xp) = 0. (5) 

In this case (1) has a parabolic fixed point p = {xp, 0) and is no longer Anosov. The main result of 
[LerlO] was to show that the diffeomorphism acquires elliptic behavior when it crosses this Anosov 
boundary. Another feature of this map is the separation of the phase space into two regions, one 
in which the dynamics is nonhyperbolic and the other in which the diffeomorphism appears to be 
nonuniformly hyperbolic. Though neither of these statements were proved in [LerlO], considerations 
in favor of these statements were presented. 

In this paper we try to use numerical methods to substantiate the following assertions about 
(1) under the assumptions (2) and (5). 

• There is an invariant, open region E C T'^ whose boundary is formed from the stable and 
unstable manifolds of two fixed points of the map, a hyperbolic saddle, h, and a parabolic 
point, p. The Lebesgue measure of E is strictly less than that of T^. 

• The channel E contains all non-hyperbolic orbits of /, and indeed has elliptic orbits for generic 

e > 0. 

• Conversely, the dynamics of fin, where H = T‘^ \ E, is nonuniformly hyperbolic; that is, the 
map is ergodic in H and has positive Lyapunov exponent. 

Of course, these statements are purely numerical observations, which should therefore be considered 
mathematically as conjectures. 

2 A Parabolic Standard Map 

Following [Lew80, LerlO], we study the dynamics of (1) using the degree-one circle map 

g{x) = x +{l + e)sm{2Trx)], (6) 

ZTT 

where g and e > 0, see Fig. 1. Note that when e = —1 and g = 0, the map (1) reduces to Arnold’s 
cat map (3). 

The map / is a diffeomorphism whenever g is smooth. Indeed 

y) = {x -y,y - g{x - y)). 

Moreover, this map is reversible, f o S = 5o/“^, with the “second” reversor of Chirikov’s map (It 
does not have the first reversor since g is not odd when g ^ 0.), 

S{x,y) = {x -y,-y), (7) 
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Figure 1: Graph of the force g (6), for e = 0.5 and /r = g,p{e) ~ 0.27696 from (9), with the parabolic point 
Xp « 0.13386 and saddle xu ~ —0.27889. 


with the fixed set Fix(5) = {(s, 0) : s G §}. Note that since S is an involution, the map 
f o S{x,y) = S o y) = {x-2y + g{x - y), -y + g{x - y)) 

is also a reversor, with the fixed set 

Fix(/ o S) = {s+ y{s), y{s) : s G S}. 

Under the assumption (5), g{x — Xp) = 0{{x — Xp)^), the map (1) has a (symmetric) parabolic 
fixed point p = (xp,0). For the case (6) this fixed point occurs at 


xp = ^ sec-\l + ^) = ^ (l - > 

when pL is chosen to be 

pL = pip{e) = (1 + e) sin(27rxp) — 2TTXp = \/e{2 + e) — sec“^(l + e) 

= ^£ 3/2 A _ + 

3 V 20 ^ ' 


( 8 ) 


(9) 


Note that since Dg{xp) = 0, the Jacobian 

Df = 


1 + Dg{x) 1 
Dg{x) 1 


( 10 ) 


at Xp has a double eigenvalue 1 with a nontrivial Jordan block. Moreover, as was shown in [LerlO], 
the first nonzero coefficient in the nonlinear normal form is quadratic in x — Xp. 

We will fix ^ = pLp{£) using (9) and think of / as a one-parameter family f^. This family also 
has a (symmetric) hyperbolic saddle fixed point h = (x/i,0) where Xh < 0 is the negative root of 
g{x), 

Xh =(l - . (11) 


TT 


60 
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This is a hyperbolic fixed point of / since Dg{xh) > 0 whenever e > 0; indeed its multipliers are 

X± = 1± 

The pair of fixed points (xp,0) and (x/i,0) are born from a degenerate saddle at the origin for 

6 = ^ = 0 . 

Additional fixed points may exist when e is large enough so that g{x) is a nonzero integer 
for some x G [— 5 ]- For example a saddle-center bifurcation creates a pair of fixed points near 
(x, y) = (—0.2151,0) when e ~ 3.603. However, we will study the case that e is much smaller than 
this so that / has only two fixed points. 


3 Invariant Manifolds and the Channel 


The region E mentioned in §1 is an open domain that is bounded by segments of the stable and 
unstable manifolds of the parabolic and saddle fixed points p and h. In a sense, the construction 
we follow is similar that of a “DA-map” on the torus [Sma67]. In that case, the unstable manifold 
of a saddle of an Anosov map is blown-up to make a channel E in such a way that it contains 
the original unstable fixed point, and its boundaries become unstable manifolds of two saddles. 
However, since (1) is symplectic we cannot blow up a single manifold; instead, the curve to be 
blown-up corresponds to the right going halves of the unstable and stable W^{0) manifolds 

of the degenerate saddle O = (0, 0) of the almost-hyperbolic diffeomorphism /o at e = // = 0. 
Near the degenerate saddle, these two right-going curves form a cusp-shaped separatrix. Since /o is 
conjugate to an Anosov map, these manifolds are dense on the torus, and they intersect on a dense 
set of homoclinic points. Moreover, at every intersection point (except for the degenerate saddle 
itself), the manifolds cross transversely [LewSO]. As a consequence, if we thicken these two curves 
we must get a Cantor-like set on the torus. Our simulations support these conclusions, though 
the situation is more complex. To avoid these complications, we do not attempt this blow-up 
construction, but instead use a bifurcation technique. 

As noted in § 2 , when e > 0 the family has two fixed points: a hyperbolic saddle h and a 
parabolic point p. Each fixed point has smooth, stable, , and unstable, IT^, manifolds; numerical 
approximations are shown in Fig. 2. The existence of the manifolds of the saddle follows from the 
standard stable manifold theorem [HPS77]. For the parabolic point, existence follows from the 
results of Fontich [Fon99]. In contrast to the saddle, the curves IT*’“(p) are one-sided and they 
start as a cusp of the form 

IF®’“(p) = {(x, a(x — Xpf’/'^) + 0(e) : x > Xp} 


with a = Ty ^[e (2 -|- e)]^^^, see App. A. As will be discussed in §5, the manifolds of the parabolic 
point are asymptotic to the right-going branches of the manifolds of the hyperbolic point; this is 
seen numerically when the manifolds are extended, as shown in Fig. 3. 

Since the right-going unstable manifolds of p and h are born out of the degenerate manifold Wq, 
one can prove that they depend continuously on e; moreover, any finite segment of these manifolds 
detached from the fixed point can be proven to vary continuously in the C^-topology. In particular, 
if T is any segment transverse to the local unstable manifold of the degenerate saddle, then the 
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Figure 2: Stable (brown) and unstable manifolds (green) for the parabolic fixed point p = (*^,0) and for the saddle 
h = {xh,0) (blue and red, respectively) for the map (1) with (6), e = 0.1, and p = ~ 0.029558. Also shown 

(dashed curves) is a 10*^ degree polynomial approximation (21) for the parabolic manifolds. The gray region is the 
channel E. 

local manifolds W'^{h) and W'^{p) that emerge will continue to transversely intersect T when £ is 
small enough. 

Since the map has the reversor S (7), and the points h and p are symmetric, their corresponding 
stable manifolds are obtained by applying S to the unstable manifolds, and the same continuity 
and transversality conditions apply. 

The boundary of the initial part of the upper half, of the channel E can be constructed 
by following the right-going manifold from h until it first intersects the transverse segment 

T = {( 5 ,?/) : y G [0, ^)}. For small enough e, this segment is transverse W^{p) as well, so follow T 
to T n bF“(p) and then continue along the parabolic manifold back to p. Connecting p to h along 
the x-axis, closes the boundary of the initial part of E~^. A similar construction for the stable 
right-going manifolds using T' = {(x,—: x G [0, |]}, gives the lower half, E~. The union is 
shown in gray in Fig. 2. The remainder of the channel is then obtained by forward and backward 
iteration. 

One has to remark that the behavior of manifolds for h and p is considerably more complex 
than that for e = 0—especially because of the formation of folds. This is more easily seen for larger 
values of e: Fig. 4 shows the manifolds for e = 0.5. The size and number of these folds on any 
finite segments grows with e. 

These folds are caused by homoclinic intersections. To see this, choose a segment of kF“(0) 
for the almost hyperbolic diffeomorphism /q that is long enough to transversely intersect W^{0) 
at a primary homoclinic point qq. Since the perturbed manifolds of p and h are C^-close to the 
degenerate manifolds when e > 0, the unstable manifolds will transversely intersect both of the 
stable channel boundaries near qq. This gives rise to a rectangle R, a portion of the channel that 
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Figure 3: Numerically computed stable and unstable manifolds for the saddle (xh,0) (left) and the parabolic point 
( 2 ;p, 0 ) (right) for the map (1) with (6), e = 0.1, and fi — 



Figure 4: Numerically computed stable and unstable manifolds for the saddle (blue and red) and parabolic (brown 
and green) fixed points for (1) with (6), e = 0.5, and /r = i-Lp{e). The manifolds for the saddle h are mostly obscured 
because they were plotted before those of p. 

is re-injected into itself, see Fig. 5. Two of the corners of R are homoclinic to p and h, and two are 
heteroclinic from/to p and h. Since iteration causes points to move monotonically along the stable 
manifolds, the corners on W^{h) monotonically tend to h, and the two on W^{p) monotonically 


7 













approach p. However, the interiors of the sides of the rectangle formed from the unstable manifolds 
are inside E, so they must move along the channel and tend to VH“(/i). Therefore, these segments 
lengthen along W'^{h), but have endpoints tied to W^{p). This makes them form a fold, like that 
in Fig. 5. 

For each such fold in the unstable manifolds, there is a corresponding fold in the stable manifolds 
obtained by applying the reversor (7). Suppose two such symmetric folds intersect. If one decreases 
e —>• +0, these two symmetric folds tend to the related points on the unstable and stable separatrices 
of the degenerate saddle. A consequence is that these two pieces of the stable and unstable manifolds 
of the saddle have to touch each other for some value of e > 0 on some point of the symmetry line 
y = 0. This gives a homoclinic tangency. Near tangencies can be seen in Fig. 4. This story requires 
a proof, of course, but shows one mechanism for the creation of a tangency. 



Figure 5: Channel for e = 0.5 and ^p(0.5). Shown are segments of the stable and unstable manifolds of h and p as 
in Fig. 2. The unstable channel, is shaded gray, and a fundamental domain, D, is blue. The curves show the first 
re-injection of the unstable manifolds into the channel and subsequent formation of a fold. The unstable manifolds 
are clipped when they first cross the stable boundary of the channel E. The transversal crossing creates a rectangle 
R (dark gray) whose corner points are heteroclinic from/to h or p. A primary homoclinic point q is labeled. Nine 
images of R are shown. Upon iteration the unstable sides of the rectangle begin to fold due to their accumulation 
onto lF“(/i). 


4 Elliptic Dynamics 

It is known that the quadratic homoclinic tangencies (as can be seen in Fig. 4) and the resulting 
Newhouse phenomena [New77, Dua08, GTS07] create elliptic orbits. However, there are no visible 











elliptic orbits in the numerical computations when e is small. As is seen in Fig. 6, when £ is not too 
large the map appears numerically to be ergodic in the sense that a single orbit lands in every pixel 
in a computer generated image. The first visible loss of ergodicity on this scale is due to a period- 
five, saddle-center bifurcation at e ~ 1.159, which creates a chain of five elliptic islands inside the 
channel. These islands undergo the usual sequence of area-preserving resonant bifurcations, before 
being destroyed by period-doubling near e = 1.274. 



Figure 6: Numerically computed stable (blue) and unstable (red) manifolds the hyperbolic point h of (1) with (6), 
£ = 0.37970035, and = fJ-pis)- The background (yellow) shows 10^ iterates of (0,0); these cover every pixel in the 
950 X 950 pixel image. Nevertheless, there is an elliptic period-12 orbit in the channel; it is covered by the black 
boxes (the boxes are much larger than the islands!). An enlargement of one of these islands is in Fig. 7. 

For smaller e, the islands are too small—in size and/or interval of stability—to be easily observed 
numerically. Indeed, for 0.01 < e < 1.15 in steps of Ae = 0.01, we observe no islands in the phase 
space whose size is larger than 10“^. 

Nevertheless, such orbits can be found using the reversibility of the map and thus exploiting its 
geometric properties. To this end, recall an assertion from the theory of reversible maps [LR98]: 

Theorem 1 (Devaney, 1976). Suppose that 0(z} is a symmetric periodic orbit of a reversible 
diffeomorphism f with a reversing involution S having a smooth submanifold of fixed points Fix(S'). 
Then: 


• 0{z) has period 2p if and only if there exists ( G 0{z) such that ( G Fix(5) n /^(Fix(5)) or 

CGFix(/oF)n/P(Fix(/o5)); 

• 0(z) has period 2p+l if and only if there exists ( G 0{z) such that ( G Fix(S')n/^(Fix(/ o S')). 

Consequently, in order to find a symmetric periodic orbit one should search among the inter¬ 
sections of the images of Fix(S) and Fix(/ o S). In the two-dimensional case, if the intersection of 
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these two sets is transverse then this orbit is a saddle or an elliptic point. If this intersection point 
is a quadratic tangency, then the orbit is parabolic, and generically corresponds to a saddle-center 
bifurcation upon varying a parameter. The proof of this statement is in the App. B. 

This indeed is confirmed by calculation.^ For example, using the reversor (7), we found a pair 
of period-12 orbits born on Fix(5) n /®(Fix(S)) in a saddle-center bifurcation at e ss 0.37970034. 
For slightly larger e, one of these orbits is elliptic, and is surrounded by an island, see Fig. 7. As is 
typical, this island is also surrounded by other elliptic points—a prominent period-12 x 7 orbit is 
visible—and this orbit also surrounded by more elliptic periodic orbits, for example one of period 
12 X 7 X 6 = 504. This island has stable orbits only for a parameter window of width Ae ~ 2(10)“'^. 

We conjecture that there are similar, higher-period, elliptic orbits for arbitrarily small, positive 
e. However, even a numerical investigation of this seems very difficult. 



Figure 7: Island around a period-12 elliptic point for (1) with (6) for e = 0.37970035 and /r = fJ,p{£)- The bounds 
of this figure are [-0.13527,-0.13505] x [-5.5(10)"®, 5.5(10)"®]. 


5 Asymptotic Behavior of the Channel 

Let F : —>■ be a lift of /, chosen so that its fixed points, which we will call h and p, lie in 

the fundamental square [—^, \) x [— j, g)- The lifted unstable manifolds, W^(h) and VF“(p), see 
Fig. 8, are observed to be asymptotic to a ray with fixed slope, and indeed to be asymptotic to 
one another. Of course, reversibility implies that the same properties hold for the stable manifolds 
(with another slope). 

To explain this, we first parameterize the manifold W^{h) in the following standard way. Take 
some point {xh,yh) = Ch & W'^{h), and let 

£ = iT“(a,F(a)) 

^the authors thank A. Kazakov for his supporting numerics made upon our request. 
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be the fundamental segment of unstable manifold connecting to -F(C/i)- By selecting Ch close 
enough to h it is always possible to have I belong to the unit square. We parameterize i = {C/i('S) : 
s G [ 0 , 1 ]} so that C/i(0) = Qh and C/i(l) = F{Ch)- This parameterization can be extended to M using 
iteration so that s + 1 corresponds to the image. 


a(s + i) = F(a(s)), 

and s — 1 to the preimage, etc. Since the iterates of I cover W^{h), we get a full parameterization: 

W^Ch) = {a(5) = {xh{s),yh{s)) : s G M,}. (12) 

Note that Ch{s) —)• /i as s —>■ —oo. 

Computing the manifolds up to x = 1000 and for values 0.01 < e < 1.0, we observe that 


yhjs) 

Xh{s) 


A-l 


+ 0 


1 


Xhis)J ’ 


(13) 


where (f) = ^(1 + \/5) is the golden mean. The same considerations apply to the unstable manifold 
of the parabolic point W'^(p). To verify these observations, we first recall that the map (1) is 
semi-conjugate to the linear, Anosov map a, (3), i.e., that there is a continuous, onto map k that is 
homotopic to the identity such that (4) holds [LerlO]. In particular, this implies that / induces a 
map f* = a* on the fundamental group of the torus that is hyperbolic. The Anosov map a has 
a unique fixed point at the origin, a saddle with eigenvalues cjP' and The right-going unstable 
manifold of the origin, is the projection of the right-going unstable eigenvector of the matrix A 
onto the torus. This has slope (f )~^—and of course this slope is precisely the slope that we observe 
in (13). 

We will prove the following result. 

Theorem 2. Let F be a lift of (1) to the plane and assume that the force g is of the form (6) with 
e small enough. Then the right-going unstable manifolds of the saddle h = {xh,yh) o,nd parabolic 
p = {xp,yp) fixed points of F both tend to oo and have a limiting slope 4>~^, the inverse of the 
golden mean. 


As a start to the proof of Th. 2, we first note that the semi-conjugacy collapses the stable 
and unstable manifolds of both p and h onto the corresponding manifolds of the fixed point of the 
Anosov map: 

Lemma 3. The semi-conjugacy k transforms W'^{h), the upper boundary of the channel E for f, 
onto the unstable manifold 7 = W'^{0) of O = (0,0) under the Anosov map a. The same is true 
for W'^{p). Similarly both W^{h) and W^{p) are mapped by k onto the stable manifold of O under 
a. 


Proof. The map / has exactly two fixed points: h, and p, and the map a has the unique fixed point, 
O. Using (4), the fe-image of a fixed point for / must be a fixed point for a, and hence k{p) = 
k{h) = O. By definition, the backward orbit of any point C £ IU“(/i) tends to h: Ct = f^iC) h as 
t —>■ — 00 . Under the semi-conjugacy, one has k{Ct) = a* o k{C), thus the fc-image of the backward 
orbit of C is the backward orbit of the point k{(). Since k is continuous, k{Ct) —^ k{h) = O as 
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Figure 8: Stable (blue) and unstable (red) manifolds of the hyperbolic fixed point, and stable (brown) and unstable 
(green) manifolds of the parabolic fixed point for e = 0.6 for the lift F. The gray region, labeled F"*", is a portion of the 
upper half of the lifted channel. Also shown are the strips constructed from the lifts of manifold segments q) 

to a homoclinic point q. The image F{q + m) lies on the stable segment lT“(h + m', g + m!) where m' = (3, 2). 

t —>■ —oo, thus k{C) G 7 - Thus kiW'^{h)) is a subset of 7 . To prove that this image is onto 7 we 
need to find a fundamental segment in W^{h) whose image is a fundamental segment in 7 . One 
way to do this is to recall that for e = 0 the map feo is a homeomorphism. Therefore since k^ 
depends continuously on e, when e is small enough, the /t-images of the distinct points C/i, feiCh) 
are distinct. Thus the image of the fundamental segment in lT“(/i) covers a fundamental segment 
of 7 . 

By a similar argument, the same results hold for the stable manifold. □ 

The main tool we will use in the proof of Th. 2 is a theorem stated by Weil in 1935 at the 
Moscow Topological Conference [Wei36], but proved much later by Markley [Mar69]. Let L be a 
continuous, semi-infinite, simple (without self-intersections) curve on the torus = "E? jJ? and 
L = {{x{t),y{t)) : t £ M+} be its parameterized lift to the plane. Let p{x,y) denote the standard 
Euclidean distance from (x, y) to the origin. Then the following theorem holds. 

Theorem 4 (Weil (1935)). If p{x{t),y{t)) tends to infinity as t ^ 00 , then L has an asymptotic 
direction; that is, either there exists a slope m* £ M such that lim.t^aoy{'t)/x{t) = m*, or, if this 
ratio is unbounded, then limt^oo x{t)/y{t) = 0. 

In other words, if the lifted curve {x{t),y{t)), t £ [0,oo), has no finite accumulation points (i.e.. 
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there is no sequence tn ^ oo such that lim^^oo p{x{tn),y{tn)) = {x*,y*)), then it has an asymptotic 
direction. 

Given these results, we now proceed to prove our theorem: 

Proof of Th. 2. Choose lifts F and K of the maps / and k to the plane such that the fixed points 
h and p of F lie in the unit square [—|, x [—and that K{h) = K{p) = O. We are to prove 
that the unstable manifolds W^^h) and Vh“(p) of F go to infinity and have an asymptotic direction 
equal to (f ~^—the same as that of the unstable direction of the linear map A. For the first part, 
we will apply Th. 4, and thus we only need to verify that both of the lifted curves have no hnite 
accumulation points. We shall prove this for W^{h), since the proof for W^{p) is the same. 

Begin by choosing a primary, transverse homoclinic point q G VF“(h) n W^{K)f^ and consider 
the segments g) and W^{h,q), with orientations from h to q. We claim that it is possible 

to choose q so that W^{h,q) lies in the interior of the unit square. Consequently, the forward 
images, f^{q), remain in a neighborhood of h, moving monotonically from q to h along the local 
segment W^{h,q). Note that the tangent vectors at q to W'^{h) and W^{h) form a frame. Since / 
is symplectic and thus orientation-preserving, the orientation of this frame is preserved under Df. 

The existence of such a, q = q^ follows from the closeness of the manifolds of /g to those of 
/o, apart from an 0{y/e)- neighborhood of O, and the existence of a transverse homoclinic point on 
the right-going unstable manifold of VF“(0) for /q. Indeed [LewSO] showed that the intersections 
of the stable and unstable manifolds of /o are transverse everywhere except at O. Choose one 
such primary intersection, r G VF“(0) C W^{0). Since /o(r’) O as t ^ oo, there is an image 
qo = f^{r), such its forward images lie in a ball of radius, say, ^ of O. This homoclinic point qo 
is of course still primary and transverse. Now we take e small enough in order that: (1) qo is not 
in an 0(y^)-neighborhood of O; (2) the intersection point qo continues to a point q^ still in the 
unit square; and (3) the intersection at qg remains transverse. This can be done since since the 
manifolds of fg are close to those of /q. 

The loop C = W^{h,q) U W^{h,q) is a simple (since q is primary) closed curve on the torus. 
This curve is not homotopic to zero and has some nontrivial representation {mi, m 2 ) G in the 
fundamental group of the torus. This follows from that fact that C is homotopic to the loop made 
up from pieces of kF“(0) and W^{0) of /o, and hence it is homeomorphic to the related loop of 
the Anosov map a, which is not homotopic to zero. 

A lift of C to the covering plane unwinds to an infinite curve that tends to infinity with rational 
slope m 2 lm\. The collection of all lifts of C cut the plane into infinite number of disjoint strips, 
recall Fig. 8. Consider the segment U = IF“(/i, q + m) that belongs to the upper boundary of one of 
these strips, say Hi. The image FifJ) = W'^{h, F{q + m)) expands and, by orientation preservation, 
enters the interior of IIi. The second endpoint, F{q m) = F{q) + m', where m' = Am, lies in the 
interior of the segment W^{h + m',q + m'), and is not on the upper boundary of Hi, since m' is 
not parallel to m. 

Two cases are possible. The hrst occurs when FifA) —U intersects the boundary of Hi only at 
F{q + m), i.e., it intersects no other lift of W^{h, q). The implication is that at the next iteration, 
F‘^{1/() will cross the strip below IIi, due to preservation of orientation, etc. In this case, the right¬ 
going manifold W^{h) goes to infinity and cannot have accumulation points in finite part of the 

^Recall that a homoclinic point g of a plane diffeomorphism with a saddle fixed point h is “primary,” if the closed 
loop W (h, q)yjW {h, q) has no self-intersection points. 
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plane. Similar considerations were used in [Gri77]. However, it may be the case that FiJA) —U 
intersects some additional lifts of the segment q) that belong to the boundary of Hi. In this 

second case, F{hl) has to leave and then return to Hi since its extreme point F{q-\- m) still exits 
Hi through W^{h + m', q + m'). This leads to a potential problem exemplified by the folds shown 
in Fig. 5: there could be loops homotopic to zero made up from pieces of W'^{h) and W^{h). 

Nevertheless, as we show in App. C, using specific properties of the map (1), all points of W'^{h) 
tend to infinity. To apply the considerations of App. C, we need to choose a fundamental segment 
£ G W^{h) in the first quadrant, such that for all {x,y) G £, 

x + (p-^y > ^(l + 2e). 

To this end, it is enough to verify that a segment io exists for /o (for e = /r = 0), a fact that 
is easily numerically verified. Then, satisfies the requirement for for small enough e, since 
the unstable manifolds of both h and p are C^-close to those of /o on compact sets away from an 
0(\/e)-neighborhood of O (in fact, C^-closeness of manifolds is sufficient). 

Therefore we have shown that there are no accumulation points, and Th. 4 applies, and thus 
W'^{h) has a limiting slope. The same argument applies to W'^{p), since it too is mapped onto T 
by K. 

Finally we need to prove that the limiting slope is indeed equal to </>~^. To that end we use (1) 
with the assumption that g has degree one: 


g{x) = x + g{x), 


(14) 


where g{x) is a continuous, periodic function. If {xfi{s),yhis)) = Ch{s), then under the map 
F{x, y) = Qh{s + 1). Define the slope of the chord from the fixed point to Qh{s) by 

, X Vh{s) - yh 

m{s) = - T-\ -• 

Xh{s) - Xh 


Subtracting the 
algebra. 


fixed point from both sides of (1), and computing the slope gives, after some 


m{s + 1) 


1 + m(s) 


g{xhis))-g{xh) 

Xh(s)-Xh 


2 + m{s) + 

W I Xh{s)-Xh 


(15) 


Now, according to Th. 4, the slope m{s) has a limit, m*. Moreover, since Xh{s) is unbounded and 
g{x) is periodic, 

gjxhjs)) - gjxh) _ 

™ Xh{s)-Xh 

Thus after taking the limit on both sides of (15) we come to 


m 


* 


1 + m* 

2 + m* ’ 


which implies, since m{s) > 0, that m* = (j) ^ 


□ 
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Not only do W'^{h) and W^{p) for the lift F have the same limiting slope, as implied by Th. 2, 
but they converge to each other. This can be seen numerically by choosing the first parameter value 
for which each curve crosses a particular abscissa value ^ = Xp{sp) = Xh{sh)- Let T]h{C) = Vhish)-, and 
^p(0 = Vpi^p) denote the corresponding ordinates. We observe (again computing up to ^ = 1000) 
that the vertical distance between these curves decreases algebraically as 

VhiC) - VpiC) = O {C^) , (16) 

see Fig. 9. A proof of this result appears to be nontrivial. Indeed, this decay is not uniform—since 
we define ^ to be the first horizontal crossing, the formation of folds causes the vertical distance to 
exhibit jumps. More generally, close approaches to the stable channel cause oscillations; the first 
place this occurs is near (^ = 2.5. In this case the vertical distance decreases monotonically up to 
^ = 2, it subsequently increases as the unstable channel crosses the initial segment of the stable 
channel, reaching a local maximum near {xh,yh) ~ (2.59,1.64). The next two local maxima occur 
at (5.39,3.38) and (7.23,4.57), again correlated with crossing the stable channel. As can be seen in 
Fig. 9, these local maxima occur at approximately the same values of ^ for any value of e. Indeed, 
the function [r]hiO ~ appears to be quasiperiodic, with two dominant periods = 2.61 

and 6.90, again independent of e. 



Figure 9: Vertical distance between the unstable manifolds of the saddle and parabolic fixed points for the lift of 
(1) up to ^ = 200 for three values of e. Also shown (dashed curves) are the graphs a/^, with the values of a as shown. 


6 Channel Area 

Our goal in this section is to compute the area contained in the channel E “between” the invariant 
manifolds of the parabolic and hyperbolic points, and to show that, when £ is small, the total area 
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of the channel is less than one, implying that the dynamics of / is partitioned into two invariant 
regions of nonzero measure. As noted above, and in particular in Fig. 6, it is numerically infeasible 
to simply iterate chosen initial conditions in the channel, since these numerical trajectories rapidly 
fill every pixel in the image. Thus, instead, we will use numerical computations of the stable and 
unstable manifolds that form the boundaries of E. 

For the lift F, the region “between” the curves W‘^{h) and W'^{p) corresponds to the upper, 
unstable half, E~^, of the lifted channel, e.g., the gray region in Fig. 8. Its area can be easily 
computed up to some finite extension ^ on the plane. Let E~^{C) C denote the region with 
boundary 

dE+{C) = {ix,0):xh<x<Xp} + W^{p,Cp) + V^-W^{hXh) (17) 

where Ch,p = points on the unstable manifolds, and 

= {(C, y) ■ VpiO <y< r]hiO} ( 18 ) 

is the connecting vertical segment. Let 

= Area(.E+(0) 

denote the area of the unstable channel up to This can be most easily computed using the 
relation between action and area, see App. D. The results are shown in Fig. 10 for several values 
of the cut-off as a function of e. 

Since, by (8) and (11), Xp — Xh = and the slope of the unstable eigenvector of (10) at Xh 

is 0(e^/^), it can be seen that A'’'(^) = More precisely, this follows from the Hamiltonian 

normal form valid near e = p = 0. This asymptotics is supported by the calculations shown in 
Fig. 10. 

Note that A+(^) must be unbounded as ^ oo: this is a consequence of area preservation. 
Indeed, the unstable channel can be generated by iteration of a “fundamental domain.” For each 
point C = (?)^)) on a branch of an unstable manifold, the segment i = iy“(C, T(())) generates the 
entire manifold under iteration. A fundamental domain, D, for the unstable half of the channel 
corresponds to the region with boundary 

dD = W^iCp, FiCp)) + F{V^) - lT“(a, F{Ch)) - v^, 

where ilh,p) are points on the respective manifolds, recall Fig. 5. Note that the image 

F(V^) of the vertical segment (18) is a line segment with unit slope since (1) has constant twist. 
The channel E~^ is generated by the images of D, and thus, with each iteration of the map, its area 
grows by the area of the fundamental domain. 

Since A+(^) —)■ oo as ^ oo, the vertical distance between the manifolds, (16), cannot decrease 
more rapidly than confirming the decay observed in Fig. 9. Given this rate of convergence, it 
is clear that A~^ must increase logarithmically with the intercept and this is confirmed by the 
computations in Fig. 11. The oscillations seen in this figure correspond to those seen in the vertical 
distance in Fig. 9. 

So, how can we conjecture that the projection E of the channel E onto the torus has finite area? 
This must happen by the creation of heteroclinic orbits that lead to re-injection of the channel, as 
we noted in §3, and showed in Fig. 5. The region R in this figure, and all of its forward images. 
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Figure 10: Area of the lift of the upper, unstable half of the channel for the lift F of (1) as a function of e for six 
values of the cut-off 



Figure 11; Area of the upper channel on the plane as a function of the horizontal extent, ^ of the channel for 
e = 0.5 (blue, left axis) and e = 0.1 (green, right axis). Fits (black) are to the log functions shown using data up to 
i = 1000 . 
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are inside the channel, and their areas are be deleted from the lifted area of the channel upon 
projection. 

To compute the area oi E = E~^UE~ on accounting for the overlap of the channel with itself, 
we resort to an image-based calculation. To start, the region E~^ is computed up to an extension 
^ as before. This region is projected into T^, and discretized into an N x N pixel image—a pixel 
is deemed to be occupied if there is a point on the manifolds or W'^{p) that lands in the 

pixel. The region is computed by filling the pixels vertically from r]p to rjh. We fill in the channel 
sequentially, increasing the cut-off in this case the folds, which are in the interior of the channel, 
do not cause a problem with the hlling algorithm. Finally the full channel is computed by applying 
the reflection S to the pixels in the unstable channel, giving the image E]\fxN- An example, for 
e = 0.5, is shown in Fig. 12. 

We observe that as the number of pixels, N, grows, the computed channel area monotonically 
decreases, and that the error is proportional to iV“^, see Table 1. We can use this to extrapolate 
to get an estimate of the area to an absolute error less than 10““^, the column labeled 67V in the 
table. A final extrapolation to remove errors proportional to N~‘^, the column cat, reduces the 
error estimate slightly. Thus we estimate that 

Ao.i(lOO) = 0.03679 ± 2(10)“^ 

Note that the area of the upper lifted channel (computed using the action) is A^j^(lOO) = 0.018563, 
which when doubled gives a total channel area of 0.037126. Thus the fraction of area excluded due 
to overlap is about 0.9%. 


N 

ON 

bw 

CN 

500 

0.486252 



1000 

0.289216 

0.092180 


2000 

0.167454 

0.045692 

0.030196 

4000 

0.103127 

0.038800 

0.036503 

8000 

0.070165 

0.037203 

0.036671 

16000 

0.053519 

0.036873 

0.036763 

32000 

0.045162 

0.036805 

0.036782 

64000 

0.040973 

0.036784 

0.036777 

128000 

0.038890 

0.036807 

0.036815 


Table 1: Area of the discretized channel, otv, for e = 0.1 and ^ = 100 as a function of the number 
of pixels, N X N, m. the image. The extrapolation, 67V = 2a7v — aN/ 2 i removes errors 0{N~^) and 
the second, removing errors is ctv = — ^7v/2)- 

After this extrapolation, we vary the cut-off, to attempt to estimate = lim^_^oo Area(£'£(^)). 
The results, again for e = 0.1 are shown in Table 2. After the second extrapolation, we estimate 
that the true area of the channel is 


Ao.i = 0.03990 ±5(10)“^ 

Using these ideas, we compute the area as a function of e for three values of the channel cut-off, 
see Fig. 13. For e < 0.01, the results have not converged: they depend on | significantly. It 
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Figure 12: Discretized channel Fjvxiv for e = 0.5 with ^ = 100 and N = 4000. The channel (black region) intersects 
7, 040, 240 pixels, or 44.00% of the area. For ^ = 100, rjh — rjp = 1.47(10)“^, so that the channel height is less then one 
pixel. Note that the folds in the manifolds are hidden since they occur in the interior of the channel, recall Fig. 4. 


is strange that for these values, x4(200) > ^(300); this is due to error in the extrapolations for 
bN —none of the computed values a^i^) have this contradictory property. The error bars in the 
fignre are estimated by 1 632000 ~ 6i6ooo|- When s > 0.01, the area seems to have converged with 
^ = 300. As in Fig. 10 the area again grows as A fit over the interval 0.01 < e < 1.0 gives 

A, = (1.02 ± (19) 

while a fit over the narrower interval 0.04 < e < 0.6 gives an exponent of 1.49 ± 0.04. When e 
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approaches 1, the power law predicts that 1, and, as can be seen in the figure, the area 

saturates at one. 


e 

“8000 

“16000 

“32000 

heooo 

^32000 

“32000 

30 

0.043929 

0.039000 

0.036549 

0.034071 

0.034098 

0.034107 

70 

0.059595 

0.047852 

0.041964 

0.036109 

0.036076 

0.036065 

no 

0.073580 

0.055382 

0.046242 

0.037184 

0.037102 

0.037075 

150 

0.087464 

0.062759 

0.050329 

0.038054 

0.037899 

0.037847 

190 

0.102100 

0.070515 

0.054570 

0.038930 

0.038625 

0.038523 

230 

0.116212 

0.077799 

0.058356 

0.039386 

0.038913 

0.038755 

270 

0.128380 

0.084259 

0.061867 

0.040138 

0.039475 

0.039254 

310 

0.141428 

0.091035 

0.065402 

0.040642 

0.039769 

0.039478 

350 

0.154256 

0.097727 

0.068905 

0.041198 

0.040083 

0.039711 

390 

0.167162 

0.104657 

0.072596 

0.042152 

0.040535 

0.039996 

430 

0.180157 

0.112024 

0.076516 

0.043891 

0.041008 

0.040047 

470 

0.192669 

0.119685 

0.080630 

0.046701 

0.041575 

0.039866 

510 

0.205416 

0.126197 

0.083939 

0.046978 

0.041681 

0.039915 


Table 2: Channel area for e = 0.1 as a function of the cut-off Columns labeled ajv are the 
computed areas for N x N pixels. The final three columns show extrapolations bi\f = 2a^ — ajv/ 2 ) 
and CN = \{^bN - bN/ 2 )- 


7 Lyapunov Exponents 

The Lyapunov exponents of the family /g appear, by the standard computation, to be positive. We 
compute the finite-time exponent 

K{x,y,T) = ^ In \\DfJ{x,y)vo\\, (20) 

for an initial condition (x, y) with the vertical initial deviation vector vq over a time T. The 
results shown in Fig. 14 give the mean exponent for 400 initial conditions with T = 10^ (for these 
parameters standard deviation of the distribution of exponents is smaller than 0.005). Note that 
(Ao) ~ 0.902177 < ln(/)^ Ri 0.9642, the exponent of the Anosov map (3) (i.e., e = —1 and /r = 0 ). 
The exponent decreases monotonically from its value at e = 0 until e = 1.55, when it begins to 
increase (though not monotonically), reaching (A) ~ 1.1 at e = 5. 

To estimate the exponent separately for orbits that lie in the channel E and orbits that lie in its 
exterior, H = we use the N x A-pixel approximation of the channel, Ej^xn{0^ recall Fig. 12. 

Initial conditions for the exterior computation are chosen in each pixel of n{^)=T^\En>cn{^), 
and each is iterated only over the time that it remains in Hp^xNiO- the time T in (20) is chosen so 
that the orbit segment from 0 to T lies in H Similarly, an in-channel, Lyapunov exponent can be 
computed by choosing initial conditions in Ei^xNiO^ iterating them only as long as they remain in 
the approximate channel. The resulting finite-time Lyapunov exponents are shown, for e = 0.1, as 
a function of their initial condition in Fig. 15 using a channel cut-off of ^ = 50. The mean exponent 
for initial conditions in H is {\o.i)h = 0.8352, while (Ao.i)_e = 0.5123. 
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Figure 13: Area of the channel E using the second order extrapolation as in Table 1, for N = 32, 000, as a function 
of e for three values of the channel cut-off, ^ as shown. Error, estimated as the difference I&32000 — fiieoool, is only 
visible when e < 0.01 


The mean exponents in H and E are also shown in Fig. 14 as a function of £. Since these 
computations are for a fixed number of pixels, N = 3000, the approximation Hjstxn will vanish 
for large enough e. Since orbits leave this gridded approximation rapidly, we do not show these 
curves for £ > 1.6. Note that the exponent for initial conditions in is a monotonically decreasing 
function of e, while that for E primarily increases. It appears that the minimum of the globally 
averaged exponent (dashed curve in the figure) corresponds to the point at which the channel area 
reaches 0(1) so that the essentially all orbits are in the channel. In principle, the globally averaged 
exponent should be the weighted average of the channel and exterior results—but this is not true 
for the figure. The reason is that the computations are carried out over different time intervals. 
The latter two are averages over the shorter time during which orbit segments remain in E or in 
H. We observe that the value of (20) increases with T; the result is that both (Ao.i)e and (Aq .i)h 
are smaller than those of the global average, which used T = 10^. 

8 Conclusions 

We have provided numerical evidence for the three conjectures of §1 for a family of parabolic 
standard maps /^ : —)■ T^, (1) with force (6), that are homotopic to the Anosov map (3), but 

which have a pair of fixed points for each e > 0, one hyperbolic and one parabolic. We showed 
that the right-going stable and unstable manifolds of these fixed points bound a channel FI C T^. 
The lift E of the channel to the plane has unstable boundaries that are asymptotic to lines of slope 
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Figure 14; Lyapunov exponents for the map (1) with (6) as a function of e. The dashed curve (blue) is the average 
exponent for a grid of 20 x 20 initial conditions, each iterated 10^ steps. The dotted (red) and solid (black) curves 
show separately the mean exponents for orbits in the exterior and interior of the channel, estimated using 3000 x 3000 
pixel image with a cut-off ^ = 50. 


the slope of the unstable manifolds of the Anosov map. Since these maps are, in addition, 
reversible, the same assertion concerning the slope is valid for stable manifolds. The height of the 
lifted channel approaches zero as x~^, which is the maximal rate consistent with area-preservation. 

• We have computed the area A^{E) for the lift using the action, and on the torus using pixel- 
based computations. We show that A^{E) < 1 when e < 1. We conjecture that there is a 
transition near e = 1 where the measure of the channel reaches one. 

• We have found elliptic periodic orbits in the channel for several values of e. These are formed 
through saddle-center bifurcations near tangencies of the stable and unstable manifolds of 
the hyperbolic point, i.e., by the Newhouse mechanism. We conjecture that there are elliptic 
orbits in the channel for arbitrarily small, positive e, and that there are no elliptic orbits in 
its complement, H. 

• We have computed finite-time Lyapunov exponents for orbit segments both in the channel 
E and in its complement, H. As e —)> 0^ it appears that the former monotonically decrease, 
while the latter limit to the exponent of the almost hyperbolic map /q. This occurs even 
though a naive numerical iteration of any given initial condition appears to fill every pixel of 
a computed image. 

We hope that these results will present convincing arguments in favor of the hypothesis that a 
generic, sufficiently smooth symplectic diffeomorphism does have a positive measure invariant set 
where its Lyapunov exponent is positive and that is it non-uniformly hyperbolic on this set. This 
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Figure 15: Distribution of the finite-time Lyapunov exponents (20) for for e = 0.1. Panel (a) shows initial conditions 
in the channel and (b) in the exterior of the channel. The channel is computed up to the cut-off ^ = 50, and discretized 
onto a 3000 x 3000 grid. 


would show the drastic difference between properties of sufficiently smooth and C^-smooth sym- 
plectic diffeomorphisms where a generic case is zero Lyapunov exponent almost everywhere with 
respect to the Lebesgue measure [Boc02]. 


23 








Appendices 

A Parabolic Manifolds 

As shown by [Fon99], a map of the form (1), with a parabolic fixed point at {xp,0) such that 
g{x — Xp) = 0{{x — Xp)‘^) has a pair of stable and unstable manifolds. In the neighborhood of the 
fixed point, these can be parametrically represented as 

W'^ip) = ^ + _ _) ( 21 ) 

under the assumption that the dynamics on the manifold is parameterized by the one-dimensional 
map (T : M —>■ M, 

s I—>■ cj(s) = s -|- /32S^ -|- /3^s^ + ... 

Demanding that this set be invariant gives a set of equations that can be solved, order-by-order, 
for the coefficients {ai,/3j}. For the case (6), the result is 

irk /bvr^ 12 -|- 11/c^ 

144 ’ ^22) 

ttA; / 6tt^ 4 -|- k‘^ 

1^3 = Y’ 9^’ 

where k = \/s{2 + e). These expansions are well-defined only for k ^ 0, requiring e 7^ 0. Note that 
since > 0, has the form of a cubic cusp. 

This expansion, while useful for small s, does not give a good representation too far from 
the fixed point. For example, the degree-10 polynomial approximations are compared with the 
numerically generated manifolds of {xp, 0) in Fig. 2. 

B Creation of elliptic points from tangency of fixed point sets 

In this appendix we present a justification of the method of finding elliptic points used in §4. 
We consider only the case of an S-reversible area-preserving map, f o S = S o f~^ for which the 
involution S has a smooth line of fixed points, Fix(S'). 

Theorem 5. Suppose that f is a area-preserving diffeomorphism that is reversible w.r.t. a 
smooth involution S, and the set Fix(S) of the involution fixed points is a smooth curve. Then if 
f = Fix(5) n /^(Fix(S')) is a point of transversal intersection, it is a point on either an elliptic or 
a hyperbolic period-2p orbit, while if is a point of quadratic tangency, it is a parabolic period-2p 
orbit. 

Proof. Since ^ G Fix(S') PI /^(Fix(S')), then ^ = S{fi) and there is a point p G Fix(5) such that 
fP{p) = Consider first p = 1. Then we have f{p) = f{f{p)) = f{f) = f{S{f)) = S{f-^{C)) 
= S{p) = p. Similarly, one has /^(^) = f,. By induction, the same is true for any p G Z. Below we 
work with p = 1 to facilitate calculations. 
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According to the Bochner-Montgomery theorem [BM46] we can take two symplectic charts: V 
near rj with coordinates (x, y) and lA near ^ with coordinates {u, v) such that in V the involution 
S becomes S{x,y) = {x,—y), and similarly in lA it becomes S{u,v) = {u,—v). Moreover, /|v = 
/i : V ^ lA is written as follows (we assume with no loss of generality that ^ and rj have zero 
coordinates in the related charts) 

J \y)^\Giix,y) 


where A is a constant matrix and Fi and Gi are 0(2). Similarly f\i( = f 2 ■ lA ^ V has the form 

- fl (A P2{u,v)\ 

Note that in both cases, du Adv = dx A dy by area preservation. 

If ^ is the point of transverse intersection of /i(Fix(iS')) and Fix(5), then two vectors (an, 021)"'" 
and (1,0)"'" are transverse, i.e., 021 / 0. In this case, when 0 < 012021 < 1, the point r/ is elliptic 
(its eigenvalues satisfy |Ai^2| = 1), while if 012021 < 0 it is an orientable saddle, and if 012021 > 1 
it is a non-orientable saddle. 

The tangency of D/i(Fix(5)) and Fix(S') at ^ implies 021 = 0 and area preservation gives 
022 = The reversibility written in both coordinate charts provides the following relations for 
direct and inverse maps /i o S' = 5 o /2 o S = S o or in coordinate form; 

. [A = -ai2\ fA -t- f ■^2(^i,-?^) ^ 

^ ’ \yj VO an J \vj Ag 2 {u,-v)J ' 


and 

«12W3:\ /F2(x,y)\ 

2 ■ 1^021 022/ Vyy V'^2(a;,?/)/ ’ 

from where we get relations: on = 622, 012 = 612, 022 = hi, hi = 0, U 2 {x,y) = Fi{x,-y), 
V 2 {x,y) = -Gi{x,-y), Ui{u,v) = F 2 {u,-v), Vi{u,v) = -G 2 {u,-v), here Ui,Vi, U 2 ,V 2 are non¬ 
linear terms of the inverse maps Denote below for brevity on = a, 012 = /3, then 

022 = a~^. 

The quadratic tangency of /i(Fix(S)) and Fix(S) at ^ implies d'^Gxjdx^ 7^ 0 at (0,0). The map 
near a 2-periodic point r\ has the form f 2 ° fi- Hence, the linear part of this map has the matrix 


2/3/a^ 


7 = 2/3/a 7^ 0. 


Let us notice that for the map near the point 7 to guarantee its fixed point be parabolic (not 
more higher degeneration) we need only to check that in the local coordinates 


xi = X + 'yy + p{x, y), yi = y + q{x, y), dxi A dyi = dx A dy 
the inequality d‘^q/dx‘^ 7^ 0 at the hxed point holds. For our case this quantity is the following 

g(O,O)=a^(O.O)-«=0(O.O). 
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From identities derived from the representation for fi and f 2 = Sof-^^oSwe get 




therefore we come to 


.2tv q 2 q ^ 


d^q, , d^Gi, , X 

0 ) - dx^ 

due to the quadratic tangency of Fix(5) and /(Fix(5)) at 


(0,0) /o 


□ 


C Orbit Bounds 


In this appendix, we obtain a sufficient condition for the forward orbit of a point under the lift F 
of the map (1) to be unbounded. This condition is used in the proof of Th. 2. 

Write the lift as 

^ xt+i \ 

where A is the matrix in (3), and g{x + 1) = g{x). The formal solution to this iteration is 


=- 4 (;;)+iiw(;). 




1=0 


(23) 


The power of the Anosov matrix (3) is easily computed in terms of the Fibonacci sequence, 


Ft+i — Ft + Ft-i, F_2 — 1, T_i — 0, 


(24) 


to obtain 


Thus (23) becomes 


A* = 


F2t F2t-1 

F2t-1 F2t-2 


= A* 


t(xo 

2/0 


t-1 




-i-i; 


j=0 


2j+l 


F 


2j 


Supposing that ||5(x)||oo = G, we can find a lower bound on the orbit as 


4-1 

xt > F2txo + F2t-iyo — G ^ .p2j+i, 

j=0 

4-1 

2/4 > F2t-lXo + F2t-2y0 — G 

j=0 

The solution to the Fibonacci difference equation (24) is 


F, = ^ 


> 


’ + 2 / , t 


(<(>‘ - 1), 
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where (p is the golden mean. Thus 


g ^ - 1+ - 1)] < 

Fy = ^ [.#“-■ 

Consequently if xq, ?/o > 0) then 

(i) 2 

Xt > [{P^\xo + p~^yo -G) -Xo- yo] , 

yt > [{p‘^^~^{xo + P~^yo -G) - Xo- yo] ■ 

Therefore, whenever 

3^0 + P~^yo > G, Xo, yo > 0, (25) 

then we have xt, yt —)• oo as t —)• oo. 

For the form (6) with y = yp(e) from (9), the sup-norm of g is 

G = ^[V £(1 + e) + 1 + e — sec ^(1 -|- e)] < ^(2e + 1). 

Thus the forward orbit of a point {xo,yo) in the positive quadrant that satisfies 

3^0 + > ^(2e + 1) (26) 



is unbounded. 

D Actions and Areas 

Areas bounded by segments of invariant manifolds of an exact, area-preserving map F : —>■ 

can be computed using the action-flux formulas of MacKay, Meiss, and Percival [MMP84, MMP87, 
Mei92]. In particular, suppose that F preserves the area form u, i.e., F*uj = u, and oj = —dv is 
an exact form. We say that F is exact, area-preserving when there exists a zero-form L : M —)■ M 
such that 

F*v -u = dL (27) 

In particular, the lift of (1) is exact symplectic with form uj = dx Ady with the Lagrangian 

L{x,y) = ^{y + g{x)f+ G{x) (28) 

where G is any anti-derivative of g. 

Suppose that z* = {x*,y*) is a hyperbolic or parabolic fixed point of F and U = W'^(z*,C) is 
the segment of the right-going unstable manifold between z* and the point C = £ W^{z*). 

Let Ct = F^iC) denote points on the orbit of (" = Coj so that Q ^ z* as t ^ —oo. 
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Consider the region TZ “below” the segment U and above x-axis, as sketched in Fig. 16. This 
region is bounded by the loop 

dn = {(x,0) : X* < X < C} + {(^,y) ■■ 0 < y < r]} - U - {{x*,y) :0 <y < y*}. 

The area of TZ is 



upon doing the trivial integrals along the straight segments of dTZ. The remaining integral of the 
one-form u along the segment lA can be done using (27), recursion, and the fact that F~^{IA) —)■ z*: 



dL + 


iF-^u) 

-1 


/ 


L(C-i) - L{z*) + 


= {L{Ct)-L{z*))^AA-{C,z*), 


f 

Jf-^{u) 


(29) 


the difference between the past actions of the two orbits. 

For the map (1), the upper half of the channel E^{^) has boundary (17). Since the fixed points 
have y* = 0, the channel area is the difference between the areas below the hyperbolic manifold 
and that below the parabolic manifold, as given by (29): 

A+{i) = Area{E+{i)) = A^-(a,h) - A^-(Cp,p). 

Areas computed using this formula for the map (1) are shown in Fig. 10. 



Figure 16; Areas below the initial segments of the right-going unstable manifolds of [xh, 0) (red), and (xp, 0) (green 
checkerboard) for e = 1.0. 


For a symmetric fixed point z*, the reversor, (7), maps S{W'^{z*)) = W^(z*). The image of the 
channel E'^ is bounded by the curve 

SidE+) = {(x,0) : X;, < X < + w^p, SiCp)) + {(^ - y,y) :-pp > y > -ph} - W^ih, 5(a)) 

Note that the reflected channel has a cut-off that is a line segment with slope minus one. Now since 
5 is area-preserving, but orientation reversing, the area of the stable channel E~ = S{E~^) is 



Thus, up to the sign, the areas are the same. 
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